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GENERALIZATION  OF  THE  SIMILARITY  EQUATIONS 
FOR  CONTINUOUS  SURFACE  BOUNDARY  LAYERS 


ABSTRACT 

Through  a  generalization  of  the  continuous  surface  problem  (i.e., 
boundary  layers  on  a  moving  belt  or  behind  a  shock  wave),  an  equation  is 
derived  which  encompasses  many  of  the  well  known  similarity  solutions  of 
boundary  layer  theory  as  well  as  some  solutions  not  previously  consider¬ 
ed.  Since  these  problems  involve  both  free  stream  and  wall  velocities, 
similarity  variables  are  introduced  which  depend  on  velocity  differences. 
The  parameter  B  =  U^/JU  -  U^)  describes  the  relative  importance  of  the 
boundary  conditions  along  with  the  usual  pressure  gradient  parameter,  B. 
This  formulation  of  the  problem  includes  the  following  special  cases: 
flat  plate  (Blasius),  accelerating  or  decelerating  flow  (Falkner-Skan) , 
boundary  layer  behind  a  shock  or  expansion  wave  (Mirels),  continuous 
surface  (Sakiadis)  and  accelerating  wall  and  free  stream  (Moore).  Two 
new  conditions  not  previously  considered  involve  reverse  flow  and  acceler¬ 
ation  and  deceleration  of  a  continuous  surface  boundary  layer.  Prelimin¬ 
ary  numerical  calculations  have  been  made  for  these  conditions. 
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I.  INTRODUCTION 


B.  C.  Sakiadis1,3’3*  published  a  number  of  papers  on  the  so-called 
continuous  solid  surface  problem  in  19bl  in  which  he  considered  the  two- 
dimensional  and  axisymmetric  cases.  In  the  simplest  two-dimensional 
situation  the  problem  is  posed  as  follows;  consider  a  smooth  flat  semi¬ 
infinite  belt  issuing  from  a  slot  irto  a  region  ot  initially  still  but 
viscous  fluid.  The  moving  belt  entrains  some  of  the  neighboring  fluid 
through  the  action  of  viscosity  and  ultimately  a  steady  state  situation 
results  in  which  a  boundary  layer  develops  near  the  belt,  beginning  at 
the  slot  and  extending  to  infinity  in  the  x  direction.  (See  Figure  1.) 

In  the  original  paper,  the  boundary  layer  equations  were  employed 
without  proof  that  they  were  applicable;  however,  the  results  of  the 
analysis  did  indicate  the  existence  of  a  boundary  layer  like  situation 
with  one  of  the  major  differences  being  that  the  y  component  of  velocity 
of  the  external  fluid  is  directed  toward  the  belt  whereas  in  the  usual 
boundary  situation,  the  y  component  of  velocity  is  directed  away  from 
the  wall.  It  was  further  shown  that  a  similarity  approach  could  be  used 
to  reduce  the  boundary  layer  equations  to  the  Blasius  equation; 

{'"  +  f  f"  =  0  f  =  f(n) 

n  =  y//2vx7TT 


with  boundary  conditions  different  from  those  of  the  Blasius  problem  of 
a  stationary  flat  plate  in  a  uniform  flow.  (Reference  1  contained  a 
factor  1/2  in  front  of  f"'  but  this  can  be  eliminated  by  a  simple  change 

of  variable.) 

Boundary  Conditions 

Continuous  Surface  Flat  Plate  (Blasius) 

f(0)  =  0  f'(0)  =  1  f(0)  =0  f'(0)  =  0 

f'(»)  =0  f'(»)  =  1 


References  are  listed  on  page  45. 
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I'ho  resulting  numerical  solutions  for  the  continuous  surface  case 
were  found  to  be  different  from  the  Hlasius  result  despite  the  fact  that 
the  equation  is  identical.  Some  of  the  differences  which  were  found  are: 

Continuous  Surface  flat  Plate  (Blasius) 

5//77IT  b.S7  4.91 

C,  'Re  0.8S0  0.&64 

f  x 


The  reason  for  the  difference  between  the  flat  plate  and  the  contin¬ 
uous  surface  solution  is  because  the  equation  is  nonlinear  and  therefore, 
the  boundary  conditions  cannot  be  transformed  without  affecting  the 
equation  itself.  The  physical  significance  of  this  mathematical  state¬ 
ment  will  be  considered  in  the  next  section.  In  the  third  section,  the 

relationship  between  the  continuous  surface  problem  of  Sakiadis  and  the 

4  5 

analysis  of  the  shock  tube  boundary  layer  of  Mirels  ’  is  discussed. 

This  leads  to  a  generalization  whereby  a  single  equation  is  derived 
which  covers  many  of  the  well  known  similarity  solutions  and  indicates 
the  possibility  of  some  new  results. 

In  the  following,  it  will  be  assumed  that  the  boundary  layer 
equations  are  applicable.  It  is  recognized,  however,  that  a  proof  is 
not  available.  Furthermore,  the  successful  numerical  integration  of 
these  equations  does  not  definitely  establish  whether  or  not  solutions, 
in  fact,  exist.  Consistent  with  boundary  layer  theory,  it  is  assumed 
that  a  characteristic  Reynolds  number  is  large  and  that  the  main  effects 
of  viscosity  are  confined  to  a  relatively  thin  region  (compared  to  the 
distance  along  the  surface)  near  the  surface. 

Consideration  will  be  limited  to  incompressible,  similarity 
problems.  The  restriction  to  similar  flows  implies  that  velocity  profiles 
when  written  in  non-dimensional  form  are  a  unique  function  of  a  non- 
dimensional  coordinate  perpendicular  to  the  surface  and  thus,  the  equations 

1  i  j 

of  motion  reduce  to  ordinary  differential  equations.  Sakiadis  ’  ’  has 
studied  this  problem  as  well  as  extending  his  work  to  include  approximate 
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integral  techniques,  investigating  the  nonsimilar  axisymmetric  and 

turbulent  situations.  Fox,  et  al  ,  have  analyzed  the  heat  transfer 

problem  on  continuous  surfaces  and  Tsou,  et  al  ,  have  measured  the  heat 

transfer  and  flow  properties  of  laminar  and  turbulent  continuous  surface 

7  8 

boundary  layers.  A  number  of  other  investigators  ’  have  also  contrib¬ 
uted  to  the  knowledge  of  continuous  surfaces  including  the  application 
to  ground  plane  simulators  for  wind  tunnels  ,  to  aircraft  flaps  and  to 
the  study  of  rapidly  accelerated  boundary  layers  .  Moore  has  also 
used  two  solutions  of  accelerating  wall  and  frecstream  flows  in  a  study 
of  unsteady  boundary  layer  separation  and  it  is  believed  that  a  more 
complete  investigation  of  the  equations  developed  here  might  also 
contribute  to  the  study  of  the  criteria  for  boundary  layer  separation. 

II.  RELATIONSHIP  BETWEEN  STEADY  AND  UNSTEADY  BOUNDARY  LAYER  EQUATIONS 

It  is  an  axiom  of  fluid  mechanics  that  dynamic  phenomena  only 
depend  on  the  relative  velocity  between  the  fluid  and  its  boundaries  and 
not  on  any  absolute  velocities,  at  least  as  long  as  all  velocities  are 
small  compared  with  the  speed  of  light.  Sakiadis  results,  however,  are 
numerically  different  from  those  obtained  using  the  same  boundary  layer 
equations  with  the  only  difference  being  a  reversal  of  the  f '  boundary 
conditions.  The  reason  for  the  difference  is,  of  course,  that  the  same 
equations  do  not  apply  when  a  uniform  velocity  is  superposed  to  change 
the  boundary  conditions.  For  example,  the  flat  plate  in  a  uniform  flow 
involves  the  steady  state  equations  of  motion.  If  a  uniform  velocity  in 
the  opposite  direction  is  imposed  to  reverse  the  boundary  conditions,  it 
is  necessary  to  use  the  unsteady  equations  of  motion. 

Consider  the  two  coordinate  systems  of  Figure  2.  Where  the  Xj,  y 
coordinates  are  fixed  in  space  and  the  plate  moves  in  the  positive  x^ 
direction  away  from  the  origin,  the  x2,  y  system,  on  the  other  hand,  is 
fixed  with  respect  to  the  plate  and  "sees"  a  free  stream  flow  in  the 
x2  direction. 
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Hie  transformation  of  any  x  coordinate  in  the  x^,  y  coordinate 
system,  into  the  Xj,  y  system  is  given  by: 


x  =  U  t  -  x 


(1) 


Differentiation  with  respect  to  time  gives  a  relationship  between 
velocities  in  the  two  systems. 


u„  =  u  -  u 


1 


(2) 


This  last  equation  must  be  considered  in  the  sense  that  U2  =  U2(x2,y) 
and  that  Uj  =  Uj  (x^y  and  t  =  constant)  where  x2  =  U^t  -  Xj.  That  is, 

equation  (2)  relates  the  steady  state  velocity  field  (subscript  2)  with 
the  unsteady  velocity  field  at  a  given  instant  in  time. 

rhe  unsteady  momentum  equation  with  consistent  boundary  layer 
assumptions  is: 


au,  an,  a2u, 

H;  *  vl  aT  '  »  5T- 


(3) 


It  will  be  assumed  that  a  solution  of  the  unsteady  momentum  equation  can 
be  obtained  in  the  form 


=  f'(n)  where  n  =  y/gUj.t) 


When  these  assumptions  are  applied  to  the  unsteady  momentum  equation, 
and  using  the  continuity  equation,  the  result  is: 


&  ii 

v  3t 


nl  fl  *  ”  8 


1L. 

3x, 


f  f  =  0 
rl  1 


(4) 
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In  order  for  equation  (4)  to  have  a  similarity  solution,  it  must  reduce 

to  an  ordinary  differential  equation  with  constant  coefficients,  and 
//  .  - 


for  the  terms  nj  fj  and  fj  fj. 

Thus 

a  =  &  IS. 

al  v  3t 

2 

implying  g  /2 

=  VUjt  +  CjtXjl 

and 

a  =  ^  a  1S- 
a2  \>  K  3X1 

2 

implying  g~/2 

=  TT  Vl  +  c: 

thus 

g2/2  =  va  t  + 

oo 

a2  X1 

and 

nl  '  y  >J2  t 

(a2xi  +  a\U<»t) 

The  differential  equation  becomes 


fl  +  {\  (Vl  +  a2fl)  ’  ° 


Consider  the  special  case  =  1  and  a2  *  _1 


f“  ♦  t\  in,  -  V  -  0 


where 


n,  = 


1  "  72  f-  (U„t  -  Xj) 


(51 

(6) 


This  is  to  be  compared  with  the  well-known  similar  solution  for  the 
flat  plate  case. 


with 


f  '  *  f,  f,  =  0 


n2  = 


•J~2  U  *7 


(71 
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liquation  (5)  can  be  transformed  into  equation  (7)  by  the  relations 

f-l  =  1  -  fj  and  n-,  =  Oj  (8)  ^ 

which  in  terms  of  dimensional  variables  are 

U-,  =  -  uj 

and 

x  =  U  t  -  x 

2  00  | 

These  are  just  the  transformation  equations  (1)  and  (2)  considered 
earlier  indicating  that  the  solution  to  equations  (5)  and  (6)  corresponds 
to  the  uniformly  translating  flat  plate  problem.  From  this  discussion 
it  may  be  concluded  that:  (a)  As  is  well  known,  the  velocity  field  for 
the  uniformly  translating  flat  plate  can  be  considered  as  the  complement 
to  the  velocity  field  for  the  fixed  flat  plate  in  a  uniform  stream  in 
the  sense  that  U.,  =  U^  -  Uj  at  a  given  instant  in  time;  (b)  The  ordinary 
differential  equation  corresponding  to  the  similarity  solution  for  the 
unsteady  velocity  field  is: 

fl  *  fl  (nl  -  V  *  0 

and  transforms  into  the  steady  case  with  =  1  -  f^'  and  n2  =  Oj  which 
are  just  the  similarity  conditions  for  the  transformation  IL  =  U  -  U, 
and  x2  =>  U^t  -  x j ;  (c)  The  fact  that  f”  =  -f*  sh^..j  that  the  steady  and 
unsteady  problems  are  dynamically  similar  (the  negative  sign  occurs 
because  the  positive  direction  of  the  x  coordinate  is  reversed  in  the 
two  cases) . 


III.  RELATIONSHIP  BETWEEN  CONTINUOUS  SURFACE 
AND  BOUNDARY  LAYER  IN  A  SHOCK  TUBE 

It  is  interesting  to  note  that  the  continuous  surface  problem  also 
has  a  counterpart  in  unsteady  flow.  If  the  coordinate  system  is  con¬ 
sidered  to  be  attached  to  the  moving  belt,  the  boundary  conditions  are 
reversed  and  the  following  equation  results: 
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Continuous  Surface 


Transformation 


Unsteady  Continuous  Surface 


f2f2  =  0 


f2(0)  =  0  ,  f2'(0)  =  1 


f2'(»)  =  0 


f|"+  (n  -  fx)fi  -  o 
f J CO)  =  0  ,  fj(0)  =  0 
f/M  =  1 


The  continuous  surface  boundary  layer  problem  is  closely  related  to 
the  boundary  layer  development  in  a  shock  tube.  When  the  coordinate 
system  is  fixed  with  respect  to  the  shock  wave  and  when  the  shock  Mach 
number  goes  to  infinity  (and  it  would  also  require  the  ratio  of  specific 
heats,  y  =  1)  so  that  the  velocity  behind  the  shock  is  zero  with  respect 
to  the  wave,  we  have  the  steady  continuous  surface  problem  as  sketched 
in  Figure  3.  The  unsteady  problem  refers  to  the  situation  where  the 
coordinate  system  is  fixed  with  respect  to  the  tube.  Physically  the 
shock  Mach  number  is  not  infinite  and  the  boundary  layer  appears  as 
shown  in  Figure  4  a  situation  which  has  been  analyzed  by  Mirels  ’  . 
Mirels,  of  course,  started  with  the  compressible  equations  which  he 
transformed  using  the  assumptions  py  =  constant  and  Pr  =  constant.  He 
has  also  considered  approximate  solutions  and  turbulent  flow,  however, 
he  apparently  did  not  consider  the  limiting  case  of  y  =  1  and  the  Ms  *  - 

IV.  GENERALIZATION  OF  BOUNDARY  LAYER  EQUATIONS 

The  following  development  is  somewhat  different  from  that  of 
references  4  and  S  in  order  to  achieve  slightly  more  generality.  Dvnamt 
similarity  suggests  that  the  usual  similarity  variables  should  be 
modified  as  follows: 


f'Cn)  = 


i.e. 


-rir  -  n  =  ifer 

(91 

♦  (Uw  -  UJ  f'(n) 

(10) 
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with  boundary  conditions  f(01  =  0,  f'(0)  =  1,  f'(»)  =  o.  For  the  sake 
of  generality  is  considered  a  function  of  x  so  that  the  boundary  layer 
equation  with  pressure  gradient  term  is  required: 


“  f  **  g 


Note  that  it  is  assumed  that  determines  the  pressure  gradient 
independent  of  the  wall  velocity  with  the  result 


///  //  it  /  9  / 

f  +  A  f  f  +  B  n  f  -  Cf  -Df  =  0 

where 

A  ■  &  4-  g  (U  -  U  ) 
v  3x  “  1  w  «' 


(11) 

(12a) 


are  the  conditions  on  g  and  U  to  insure  similarity  solutions. 

The  most  general  conclusions  which  can  b#  drawn  from  these  equations 
are  obtained  by  forming  the  following  combination: 


constant 


This  implies  that  either  (a)  Uw  is  proportional  to  U_,  (b)  U  and  U  are 
constant,  or  (c)  =  0  or  *  0.  Case  (a)  is  the  least  restrictive 

condition  and  the  other  two  cases  can  be  obtained  as  special  conditions 
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-  „  nn  I.  i]  and  x  can  be  obtained  as  follows, 
of  (a).  The  dependence  of  g  on 

First  fora: 


» -  §  ■  i  (f1 

|1  =  (A  -  §)  vx/(UK  -  UJ 


(13) 


This  equation  can  be  combined  with  the  C  equation 
dependence  of  Uw  -  U.  on  x,  provided  that  A  -  ?  * 


to  determine  the 
o,  i.e.: 


C/2 


(U  -  UJ 

V  W 


r-  C77  ■  uw  -  u„  3x  w 


so  that 


and 


„  u  =  c  xC/(2A'C)  -  C.xm 

u  “  U.  -  1 


_  2m  » 
C  =  f  ft 
^  m+1 


(14) 

(15a) 


Remembering  that  UJO^  -  U 
obtain 


)  is  a  constant,  it  is  possible  also  to 


(15b) 


B  “  if  -  x 

w  “ 


U»  r  _  4m  -  A 

0=2  U  c  rn+1  Uw  -  U. 


(15c) 


The  differential  equation  becomes 


r-AUff^ir^rnf")  -  &  cf' 


f')}  =  0  (16) 
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For  the  special  case  C  =  2A 


4 


B  = 


U  -  U 
w  “ 


c  =  :a 


(17a) 


(17b) 


D  =  4 


U  -  U 
w  " 


(17c) 


The  differential  equation  is: 

U 


f +  A  {  ff  '  + 


U  -  U 
w  “ 


nf"  -2  ( f /2  +  2  — 7T  f  ) 1  =  0  (18) 

W 


which  corresponds  to  the  limiting  acceleration  condition  when  m  The 

similarity  variable  n  =  y/g(x)  has  a  different  form,  i.e.. 


2vL 

=  u-rr 

w  ” 


(19) 


where  L  is  a  length  scale  factor  and  the  velocity  distribution  is 


U  -  U  =  (U  -  UJ„  e 


A  — 

A  L 


u»  -  u-rir  cuw  -  u-]o e 

w  00 


A  — 

A  L 


«)q  e 


A  — 

A  L 


(20a) 


(20b) 


(20c) 


where  (U  -  U  ).  is  a  velocity  scale  factor. 
W  00  u 
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It  is  usually  possible  to  eliminate  one  of  the  constants  in 
equations  (16)  or  (18)  by  a  transformation  f  =  af^  n  =  an^  however, 
if  the  greatest  generality  is  to  be  retained  A  must  change  sign  when 
U  +  U  changes  sign  although  its  magnitude  can  be  made  arbitrary  by 
the  indicated  transformation.  As  long  as  +  Uw  >  0,  the  boundary 
layer  growth  is  in  the  +  x  direction  even  though  Um  may  be  negative, 
but  when  +  Uw  <  0  the  boundary  layer  must  be  considered  to  grow  in 
the  -  x  direction  even  if  U  is  positive.  Figure  5  provides  a  schematic 
interpretation  of  the  above  concept,  if  this  concept  seems  too  artifi¬ 
cial,  one  can  always  divide  the  flow  regimes  into  two  problems  depending 

on  whether  U  +  U  is  positive  or  negative, 
w  00 

Nevertheless,  there  is  some  advantage  in  programming  for  a  computer 
to  writing  a  single  equation  encompassing  as  many  boundary  conditions  as 
possible;  thus,  the  value  of  A  will  be  assumed  as  follows: 


A  = 


U  +  U 
w  “ 


2B  +  1 
| 2B  +  l| 


(21) 


The  equation  for  case  (a)  can  now  be  summarized 


U  +  U 

w  "  (ff"  + 


f  + 


,U  +  U  | 
|  w  "I 


U»  .  2m  ,f'2  U” 

nf  "  m+r  (f  +  2  u~nr 


U  -  u 

W  1 


f')}  =  o  (22a) 


f(0)  =  0,  f  (0)  =  1,  f  (-)  =  0 


(22b) 

(22c) 

(22d) 
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(22e) 


The  special  case  (b)  where  U.  and  are  constant  is  characterized  by 
m  =  0.  The  differential  equation  reduces  to: 


U  +  U 

W  “ 


f  (f  ♦ 


u  -  u 

W  “ 


n)  =  0 


with  only  one  parameter  uyit^  -  U J  .  The  various  flows  which  are 
represented  by  this  equation  are  indicated  in  Figure  6.  The  area  of 
the  diagram  in  which  Uw  and  Uro  are  both  positive  and  U„  <  Uw  represents 
the  class  of  flows  considered  by  Mirels4*5  and  corresponds  to  the  situa¬ 
tion  of  boundary  layer  growth  behind  a  shock  wave.  There  are  two  bound¬ 
aries  to  the  region,  namely  where  there  is  no  relative  velocity  and 
where  U  =0.  The  latter  case  is  the  continuous  surface  problem  con- 
sidered”by  Sakiadis1.  Below  the  positive  Uw  axis  the  free  stream 
velocity  is  in  the  opposite  direction  to  that  of  the  wall,  i.e.,  reverse 
flow.  It  is,  of  course,  not  obvious  that  solutions  to  equation  (23) 
exist  in  this  region.  The  equation  is  singular  at  the  point  -  -Uw 
because  (U  +  U  )/|u,  +  Uj  is  discontinuous  there.  The  reverse  flow 
region  has^a  boundary  where  Uw  =  0  and  this  corresponds  to  the  flat 
plate  problem  but  in  the  unsteady  form  in  which  the  boundary  conditions 
are  reversed.  In  the  half  quadrant  where  both  Uw  and  U.  are  negative 
but  U  <  U  ,  we  have  a  slip-flow  like  regime.  It  differs  from  the  usual 
slip-flow  problem  in  that  the  boundary  condition  at  the  wall  is  a  con¬ 
stant  velocity  whereas  in  the  low  density  slip-flow  problem  the  wall 
velocity  is  assumed  proportional  to  the  local  velocity  gradient  and, 
therefore,  a  function  of  x.  The  slip-flow  situation  described  here  was 
also  analyzed  by  Mirels5  which  he  associated  with  a  weak  expansion  wave. 
The  pressure  gradient  problem  is  characterized  by  m  t  0.  The  well- 
known  Falkner-Skan  flows  are  obtained  when: 
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-1 


-1 


B  = 


U. 


U 


w 


U 

00 


iU  +  u 

w  “ 


although  the  equation  is  in  the  transformed  coordinates  appropriate  to 
the  unsteady  problem 


...  „  „  ?  . 

f'  -  if  f  ■  n  f  -  ^  (f  -  f  )}  =  o 


(24) 


with  the  boundary  condition 

f(0)  =  0,  f'(0)  =  1,  f'(»)  =  0 


The  usual  form  is  recovered  through  use  of  the  transformation  f'  =  1 


As  pointed  out,  a  number  of  the  flow  regimes  have  been  explored  in 
papers  already  in  the  boundary  layer  literature;  however,  by  presenting 
the  problem  in  the  generalized  way  used  here,  it  is  possible  to  identify 
two  regimes,  i.e.,  the  reverse  flow  and  the  Uw«  xm  cases  which  have  not 
been  considered  previously. 


V.  SOLUTIONS 

Equation  22a  has  been  programmed  for  the  U.  S.  Army  Ballistic 
Research  Laboratories'  BRLESC  computer  using  a  simple  prediction- 
corrector  integration  scheme,  A  step  size  of  .005  in  n  was  used  and 
the  integration  was  carried  out  to  an  n  of  20.  The  usual  procedure 
was  used  of  estimating  the  boundary  condition,  f"(0),  and  then  correct¬ 
ing  the  estimate  until  f'(20)  =  0  within  .00005. 

The  results  of  the  numerical  integration  provides  f  and  its 

derivatives  as  functions  of  n.  The  results  are  presented  graphically 

in  Figures  7  through  15.  The  boundary  layer  parameters  f"(0),  e(°»), 

* 

6  and  v(»)  are  given  in  Table  1. 
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« 

It  is  useful  to  indicate  the  relationship  of  the  calculated 
quantities  to  the  conventional  boundary  layer  parameters.  Thus,  we 
define  the  fol  lowing  quantit  ies: 

Shin  friction  coefficient-- 


du 

d> 


(  U  - 

v  w 


1 

x  2A 


f"  (0) 


(25) 


Momentum  thickness  — 


/OU 

U 

o 


(U  -  U)  dy  =  -g  6(») 


(26) 


Nondimens  ional  momentum  thickness  — 


Displacement 

6 


thickness 


(U 


U)  dy  =  -g  f(”)  =  -g 


.  * 

6 


(27) 


(28) 


The  integral  momentum  equation  can  be  used  to  provide  a  relation¬ 
ship  between  f*  (0)  and  9(»)  and  f(«)  +  .  The  integral  momentum  equation 
can  be  written: 

X  0  cr>  dll  0 00 

■f-  '  tc  /  0  U  (U,  -  0>  dy  .  JT  /  0 
Using  the  above  definition  for  cf,  B  and  5*  this  equation  becomes: 


t  Alternatively,  equation  (27)  can  be  integrated  using  equation  (22a) 
to  eliminate  f'“  with  identical  results. 
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+  (26  + 


1 


(29) 


_  dO 
2  ”  dx 


U  -  U 


6  ) 


U  -  U 


d  (UK  -  UJ 
dx 


The  differences  between  this  integral  equation  and  its  usual  form  are 
because  of  the  special  definitions  used  and  the  fact  that  Um/ -  11, , 
constant.  Manipulations  reduce  equation  (29)  to: 


f"(0)  =  -A 


1  +  3m 
1  +  m 


BM 


1J  -  U  m  +  1 
w 


(30) 


When  m  =  0  there  is  a  particularly  simple  relationship  between  f  (0)  and 
6  and  the  following  well  known  result  is  obtained: 


cf  ' 


£ 

x 


(31) 


The  integral  e(«°)  was  also  numerically  calculated  with  each  inte^ 
gration  of  the  basic  equation. 

In  addition  to  the  characteristics  indicated  above,  it  is  of 
interest  to  calculate  the  y  component  of  velocity,  V,  since  it  is  an 
important  indicator  of  the  behavior  of  the  flow.  It  is  obtained  by 
integration  of  the  continuity  equation: 


V  (y=0)  =  0 


(32) 


V  =  -g 


(U  -  U 

v  W  _ “ 

3x 


[n  u^-rir  +  f]  +  (1,W  -  ^  f'  -  fl 


which  can  be  reduced  to: 


■  +  i  uw  -  Um  g 


EL  =  -A  (f  + 


iTT  r^r  "  -  (33) 


For  large  values  of  n,  this  equation  goes  to  the  following  limiting 


form: 
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lim  V 

'i  — 


-A  t  fl-l  + 


(34) 


2m  U~ 

ITT]  ITT  r  nl 

W 

rhis  equation  combined  with  the  II  =  cx  satisfy  the  continuity 
equation  and  may  he  considered  the  inviscid  free  stream  velocity 
components 


< 


Vrj  --  -A  { f  f  • ) 


2m 

m  +  1 


U 


ilkiure  '  shows  the  nondimensional  shear  stress  at  the  wall,  -f"(0) 
as  a  Mjnct  <rn  of  B  =  U  /U  -  U  with  B  =  2m/m+l  as  a  parameter.  The 

'■*’  W  00 

m  0  cut vc*  show  that  - f  "(0)  decreases  as  the  parameter  B  approaches 
reverse  flow  regime  from  above  tnd  from  below.  However,  the  two 
sections  of  the  curve  are  not  symmetric  about  the  singular  point  B  « 
-1/2.  A  few  tentative  calculations  were  made  in  the  reverse  flow 
region  for  m  =  0  and  they  suggest  a  continuous  decrease  in  -f*(0), 
approaching  :ero  near  B  *  -1/2.  However,  it  is  not  clear  that  these 
calculations  represent  real  solutions  because  of  the  limited  range  of 
integration.  In  fact,  when  solutions  even  closer  to  B  >  -1/2  were 
attempted,  the  calculations  did  not  converge.  Figures  8  and  9  show  the 
x  and  y  velocity  components  (nondimensional)  for  B  =  0  and  for  B  2  0 
and  Figures  12  and  13  show  the  corresponding  velocities  for  the  B  <  -1/2 
region. 

Only  a  few  preliminary  points  were  calculated  for  the  B  4  0  case. 
Almost  all  these  calculations  were  made  for  the  B  =  .25  (m  *  1/7)  or 
6  =  -.1  (m  =  -1/4)  cases.  The  shear  stress  -f"(0)  is  higher  for  the 
B  -  .25  situation  shown  in  Figure  7  and  lower  for  the  B  =  -.1.  The 
velocity  distributions  are  shown  in  Figures  10  and  11  for  B  =  +1  and  in 
Figures  14  and  15  for  the  B  =  -1.  For  comparison  some  of  the  Falkner- 
Skan  solutions  corresponding  to  B  =  -1  and  -.199  <  B  <  2  are  shown  in 
Figures  14  and  15.  Although  they  are  not  shown  here,  the  Falkner-Skan 
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equations  in  the  adverse  pressure  gradient  case  are  non-unique;  the 
second  branch  solutions  which  arc  not  shown  also  have  a  reverse  flow 
character  with  positive  f  "(01  and  a  region  of  the  flow  near  the  wall 
moving  in  a  direction  opposite  to  that  of  the  main  stream.  The  reverse 
flows  described  by  -1  <  B  <  0  in  this  report  are  caused  by  the  wall 
moving  in  a  direction  opposite  to  the  fre.  stream  while  the  second 
branch  Falkner-Skan  solutions  have  '1^  i'.  The  second  branch  hilkner- 

Skan  are  also  a  direct  continuation  of  the  adverse  pressure  gudient 
solutions  beyond  the  "separation"  condition  corresponding  to  f  (*>)  -  0 
with  B  =  -.199.  All  the  second  branch  solutions  as  well  as  the  present 
reverse  flow  solutions  cannot  be  assumed  to  represent  physically  real 
flows.  To  be  realized,  they  must  refer  to  a  finite  body  with  effective¬ 
ly  two  leading  edges,  one  at  x  =  0  and  the  other  down  stream.  These 
solutions  perhaps  apply  in  an  intermediate  region  with  a  region  around 
the  second  leading  edge  omitted  from  consideration  just  as  the  region 
around  the  forward  leading  edge  must  be  omitted  because  the  boundary 
layer  equations  do  not  apply  there  (as  they  also  do  not  apply  at  the 
separation  point).  This  interpretation  is  quite  similar  to  the  sugges¬ 
tion  made  by  Libby  in  interpretating  che  overshoot  profile  solutions 
which  were  recently  found  for  values  of  B  <  -.199  (reference  12).  No 
attempt  has  been  made  experimentally  to  confirm  or  refute  this 
speculation,  however,  because  of  the  difficulties  in  producing  such 
flows.  The  reverse  flow  situation  with  a  continuous  belt  might  prove 
somewhat  easier  to  investigate  experimentally  and  thus  might  contribute 
to  our  understanding  of  similar  solutions. 


The  0  /  0,  non-reverse  flow  solutions  may  be  interpreted  as 
boundary  layer  solutions  behind  accelerating  shock  waves  with  the  assump. 
tion  that  compressibility  effects  can  be  eliminated  by  a  suitable  trans¬ 
formation  or  idealization  of  the  fluid  (eu  •  constant)  and  a  restriction 
on  the  pressure  gradient  condition  in  order  to  preserve  the  incompres¬ 
sible  fluid  requirement. 
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VI  .  SUMMARY 


Hv  introducing  similarly  variables  in  velocity  and  normal  coordinate 
of  the  form, 


f  '  ( n  1 


li  - 

n  rir 

w 


n  =  y/g(x) 


1 


(U  +  U  ) 
w  «■ 


the  following  general  differential  equation  can  be  derived  from  the 
incompressible  laminar  boundary  layer  equations: 


II  ♦  U  / 
< 

W  ®|  ' 


(f"f  ♦ 


u  -  D 

W  “ 


f"  n) 


rW(f/2  +  2u-4TTf'’)  =  <’ 

W  •  ' 


with  boundary  conditions  f ( O')  =  0,  f'(0)  *  1,  f'(°°)  =  0  and  Uj*  Uw«  xm. 
This  equation  encompasses  as  special  cases,  the  following  well-known 
similar  solutions: 


U„/(UW  -  UJ  8 


Blasius  (Flat  Plate)  -1 

Falkner-Skan  (Accelerating  -1 

Flow  on  a  Flat  Plate) 

Mirels  (Boundary  Layer  Behind  -1  >>  -°° 

a  Shock  or  Expansion  Wave)  .  2  <<  » 

Sakiadis  (Continuous  Surface)  0 


0 

-.199  «  2 
0 

0 


Moore  (Accelerated  Continuous  -3/2,  -10/9  <0 

Surface) 


Formulating  the  problem  in  this  way  shows  that  in  addition  to  the 
solutions  already  obtained  there  are  a  number  of  ranges  of  the  parameters 
UJ  (U  -UJ  and  8  which  have  not  previously  been  investigated.  Two  of 
these  regions  are: 
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(a)  Uw  «  xm,  U„  =  0 


(b)  -1  <  u„/uK  -  U.  <  o 


(reverse  flow) 

A  number  of  numerical  solutions  have  been  computed  and  are 
compared  with  the  well-known  solutions. 
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Table  1.  Boundary  Layer  Parameters 
12J - e  l  =f(-»)  y  fa>)  H=g*/e 

S  *  0 


Comments 


4.0 

0 

-1.718 

1.718 

.375 

-  .375 

.218 

1.36 

3.0 

0 

-1.521 

1 . 520 

.424 

-  .425 

.279 

1.55 

2.0 

0 

-1.295 

1.294 

.502 

-  .502 

.388 

1.83 

1.0 

0 

-1.020 

1.019 

,b47 

-  .647 

.635 

2.42 

0.333 

0 

-  .782 

.782 

.868 

-  .868 

1.110 

Ref.  5 

0.20 

0 

-  .725 

.725 

.950 

-  .950 

1.310 

.. 

Ref.  5 

0.00 

0.10 

0 

-  .626 

.630 

1.164 

-1.163 

1.848 

5.84 

Continui 

Surface 

0 

-  .561 

.  566 

1.682 

-1.666 

2.972 

18.51 

-0.80 

-0.90 

-1.0 

-1.333 

-2.0 

-3.0 

-4.0 


-  .287  -  .288 

-  .392  -  .393 

-  .4695  -  .470 

-  .661  -  .661 

-  .930  -  .930 

-1.226  -1.226 
-1.463  -1.463 


No 

Solutions 

1.788  -6.208  4.97 

1.411  -3.593  4.20 

1.218  -2.591  3.75  Flat  Plate 

.905  -1.369  -  Ref.  S 

.662  -  .712  2.23 

.509  -  .415  1.76 

.429  -  .293  1.50 


■1.655  1.145 

■1.282  .909 


.390  1.75 

.641  2.33 

.648  2.23 

.631  2.47 

5.98 
1.818  5.78 


No 

Solutions 
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Table  1.  Boundary  Layer  Parameters  (Continued) 


B 

8 

f"(0) 

e 

£*  =  f(~)  V  (■») 

n=5*/e 

6 

Comments 

8  *  0 

-1.0 

-  .10 

-  .319 

-  .515 

1.444 

-2.804 

4.03 

-1.0 

-  .199 

0.0 

-  .576 

2.349 

-4.078 

4.99' 

-1.0 

.25 

-  .732 

-  .396 

.945 

-2.386 

3.29 

Falkner- 

-1.0 

.50 

-  .928 

-  .350 

.804 

-2.297 

2.99 

.  Skan 
Solutions 

-1.0 

1.0 

-1.233 

-  .292 

.649 

-2.223 

2.63 

-1.0 

2.0 

-1.687 

— 

-- 

— 

2.23 

-2.0 

-  .10 

-  .758 

-1.004 

.723 

-  .720 

2.32 

-2.0 

.25 

-1.279 

-  .799 

.559 

-  .700 

2.06 

-3.0 

.25 

-1.652 

-1.059 

.436 

-  .399 

1.62 
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Figure  1.  Continuous  Surface  Boundary  Layer 


Figure  2.  Coordinates  for  a  Moving  and  Fixed  Flat  Plate 


Figure  3.  Continuous  Surface  Boundary  Layer  as  Limiting  Condition 
of  Shock  Tube  Boundary  Layer  (M  =  »,  Y  =  1) 
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Figure  4.  Shock  Tube  Boundary  Layer  (Mg  <  =0 


Figure  5.  Direction  of  Boundary  Layer  Growth 
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Figure  9.  V  Versus  n,  B s  0,  g  =  0 
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Figure  12. 


f'  Versus  n,  B  <  -l/2»  P  "  ® 
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Figure  13.  V  Versus  n,  B  <  -1/2,  B  =  0 
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< 


6.0 


Figure  14.  f'  Versus  n,  B  =  -1,  -.199  s  B  s  2 
(Falkner-Skan  Solution) 


42 


43 


Fipire  15.  V  Versus  n,  R  =  -1,  -.199  <  0  s  2  (Falkner-Skan  Solution) 
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Through  a  generalization  of  the  continuous  surface  problem  (i.e.,  boundary  layers  on 
a  moving  belt  or  behind  a  shock  wave),  an  equation  is  derived  which  encompasses  many 
of  the  well  known  similarity  solutions  of  boundary  layer  theory  as  well  as  some 
solutions  not  previously  considered.  Since  these  problems  involve  both  free  stream 
and  wall  velocities,  similarity  variables  are  introduced  which  depend  on  velocity 
differences.  The  parameter  B  *  U^/CU^  -  U^)  describes  the  relative  importance  of  the 

boundary  conditions  along  with  the  usual  pressure  gradient  parameter,  6.  This 
formulation  of  the  problem  includes  the  following  special  cases:  flat  plate  (Blasius), 
accelerating  or  decelerating  flow  (Falkner-Skan) f  boundary  layer  behind  a  shock  or 
expansion  wave  (Mirels),  continuous  surface  (Sakiadis)  and  accelerating  wall  and  free 
stream  (Moore).  Two  new  conditions  not  previously  considered  involve  reverse  flow 
and  acceleration  and  deceleration  of  a  continuous  surface  boundary  layer.  Preliminary 
numerical  calculations  have  been  made  for  these  conditions. 
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